Decoherence in a single trapped ion due to engineered reservoir 
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The decoherence in trapped ion induced by coupling the ion to the engineered reservoir is studied 
in this paper. The engineered reservoir is simulated by random variations in the trap frequency, and 
the trapped ion is treated as a two-level system driven by a far off-resonant plane wave laser field. 
The dependence of the decoherence rate on the amplitude of the superposition state is given. 
PACS number(s):42.50.Vk, 03.65.-w,03.67.Lx 
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According to quantum mechanics [1], a system can ex- 
ist in a superposition of distinct states, whereas these 
superposition state seems not to appear in the macro- 
scopic world. One possible explanation of this para- 
dox [2] is based on the fact that systems are never com- 
pletely isolated but interact with the surrounding envi- 
ronment, which contains a large number of degrees of 
freedom. The environment influences the system evolu- 
tion which continuously decoheres and transforms sys- 
tem superposition into statistical mixtures which behave 
classically [2, 3]. There are many assumptions involved in 
modeling the coupling of the system to its environment. 
For example, the nature of the coupling of a system to 
its environment is generally taken to be a linear [4] or 
a nonlinear [5] function of the position operator of the 
object. Assumptions are also made about the environ- 
ment. One example is to treat the environment as a 
reservoir of quantum oscillators, each of which interacts 
with the quantum system in question. Such a environ- 
ment is extremely difficult to control because of the lack 
of knowledge about the environment and its coupling to 
a system. A recent work [6] extends the investigations of 
the decoherence beyond the ambient reservoirs and en- 
gineer the state of the reservoir, as well as the form of 
the system-reservoir coupling. In the work, they apply 
noisy potentials to the trap electrodes to simulate a hot 
reservoir, the system of trapped ion then feels a control- 
lable reservoir, through which quantum superpositions 
are decohered into a state behaves classically. 

In this letter, we present a theoretical study on deco- 
herence of quantum superpositions of a single trapped 
ion through coupling to the engineered reservoir which 
is simulated by variating the trap frequency. We tackle 
the variations by treating them as a white noise, which 
results in decoherence of the trapped ion as you will see 
further . Before introducing a engineered reservoir, we 
consider a single ^Be^ ion confined in an rf(paul) trap 
and interacting with a plane wave laser field. The Hamil- 
tonian of this system may be written as (with fi = V) 

H^^ + V{r) + '^{\e){e\-\g){g\) 
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ion with two levels |e) and \g) in a trap V{r) = ^mw^r^, 
whereas the last two terms denote the coupling of the 
ion to the plane wave laser field. Based on this kind 
of models, nonclassical motional states such as thermal, 
Fock, coherent, squeezed and Schrodinger cat state are 
created[7,8], it opens novel possibilities to study deco- 
herence of quantum superposition[6]. For the plane wave 
excitation along the x-axis, the motional effects along 
the y and z directions are unperturbed, in this case our 
Hamiltonian eq.(l) can simplifies to a one dimensional 
model 

H^^ + V{x)+^-i\e){e\-\g){g\) 

+ ^(Ie)(.g| + \g){e\)+p,kU\e){e\ - \g){g\), (2) 
where a unitary transformation 



Ui=e 



|5)(g|+ei(--*-'=--)|e)(e| (3) 



\g){e\+H.c, 
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where the first three terms describe the free motion of the 



has also been made. The detuning S = uJeg—coL, the recoil 

shift ^ is ignored here. The last term in eq.(2) results 
from Doppler effect, which leads to the coupling between 
the nearest neighbor motional states [9] and the decoher- 
ence of a single trapped atomic/ionic qubit[10]. In the 
following, we are concentrating on the situation where 
the detuning S and the coupling strength fl^ is large, es- 
pecially, the coupling strength fi^ is much larger than 
the coupling of the system to ambient reservoir, the lat- 
ter condition is relevant to the recent engineered environ- 
ment experiment for avoiding the ambient environment 
which leads to the natural decoherence [6]. Although the 
coupling length('--^ KHz) is large, it is smaller than the 
trap frequency('--^ MHz) so that the internal spin dy- 
namics and the external motion of the ion in the trap 
occur on two different time scales. Therefore, it is useful 
to go to a rotating frame that eliminates the fourth term 
in eq.(2), which describes the Rabi oscillations between 
the two internal states. In the rotating frame, we will 
be able to understand more clearly how the last term in 
eq.(2), which couples the motional and the internal dy- 
namics, leads to decoherence when the trap frequency lo 
is modulated to be random. 

We go to the rotating frame by making a unitary trans- 
formation 



1 



(4) 



where a^: = \g){e\ + \e){g\, = |e)(e| - \g){g\ are pauli 
matrices. In the rotating frame, the system evolves ac- 
cording to 



(5) 



where is given by 



+ PxkLx{ax{t)ax + ay{t)ay + a^(t)(7^). (6) 
The time varying coefficients ax{t), cty{t), az{t) are given 

by [11] 

o s 

ax{t) = ^^(1 - cos{flet)), 



OLy{t) = ^sin{Q.et), 



Ui{t) = e-''^"*"*(|5)(g|e^(*)e'4(*)'^^e'^W'' 



+ |e)(e|e^^*^e 



(9) 



where 



A{t) = '-^{e'^' - l),B{t) = -A*{t), 



We consider a initial state of the form in the rotating 
frame 



IV^(O)) = Cg\g) \ag) + Ce|e) O \ae), 



(10) 



where \ai), {i = g, e) denotes a coherent state, Cg and Ce 
are constants and satisfying \cg\^ + \ce\^ = 1. This kind 
of states may be created by during Raman transitions [7]. 
Eqs(3),(4) and (9) together govern the evolution of the 
system. With these equations, we can analytically evolve 
the initial state (10) to obtain 



J2 ^2 



where Qe = + what follows, we make two 

simplifications in order to extract out the dominant be- 
haviors of the system. Noticing the coefficients ai{t),i = 
X, y, z oscillate rapidly, we expect the system in the rotat- 
ing frame to evolve on a much slower time scale than the 
period 27r/f2e- In this sense we can simplify the Hamil- 
tonian Hi{t) given in eq.(6) by taking a average values 
of the coefficients ai(t),i = x,y,z, this is equivalent to 
coarse graining eq.(6). The coefficients ai{t) then become 
time independent and reduced to dx = oiy = 0, 

otz = Furthermore, we make assumption that the 
system is being driven far-off-resonance, i.e, S » CIl- 
We therefore have 



Hi 



(7) 



we note Px 



^(a^ — a) with a(af) the annihilation 



(creation) operator for motional state \n) which satisfies 

(lm + = + 5)1^)1 therefore the Hamilto- 

nian(7) can be rewritten as 



Hi = ujo'a + ig(Tz{a) — a). 



(8) 



where g = ^^fpkLx- Hi couples nearest neighbors of 
motional state \n) with a same internal level \g) or |e). 
As will be clear further, the second term in Hi leads to 
decoherence while the engineered reservoir applied. To 
show this, wc first of all give the time evolution operator 
in the rotating frame. 



|V;(t)) = ei--*--«*«*(ai(i)c,(i)|a,-(0) 
+ a2{t)ce{t)\a-m®\g) 



+ e 
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'^\a2{t)cg{t)\a+{t)) 



+ al{t)ce{t)\at{t))) |e). 



(11) 



where 



ai{t) = cos—t - 



— sin—t,a2{t) = -i—sin—t, 



Ce{t) 



c„e 



f{t)-^\A(t)\' 



ofit)-k\A(t)[^ 



(^t{t)=ag + A{t)±-kLx 



2mu) 



-iujt 



-iujt 



a, {t)=a,-A{t)±-ki.x^^^. 

Eq.(ll) is the main result of our study, with which wc can 
explain the essential properties of the system. As eq.(ll) 
shows, the initial motional state \ag) and |ae) are dis- 
placed to be Ictg (i)) and |a^(i)), respectively. Their dis- 
placements depend on the internal state of trapped ion. 
We are now interested in the coefficient of the off-diagonal 
element of the density operator p{t) = |V'(i))('i/'(OI '^^ 
ionic internal space, its module which represents qualita- 
tively decoherence of the system is 



2 



R{oj,t) = Mo4at(0a2(t)|cg|'(a;(t)|a+(i))e-^"<(*K-*(*)* 
+ af{t)c;ce{a-{t)\at{t))e-''^'^tit)a-it)t 



+ |a2Wpc:c,(a-(0K(i))e-^-«="W<*«* 
+ at{t)aUt)\ce\'{a-it)\atit)) 

■ ^-ioja- (t)a+* (t)t |g-4^ (l-cosa;t) _ 



(12) 



Where Mod{...} denotes the module of the term in the 

brace. Wc take a special case ag = = 2 with fixed 
(J = 27r X 11. 3M Hz as an example to illustrate R{uj,t) 
vs. time in figure 1. As figure 1 shows, R{u!,t) is a pe- 
riodic function of time t, which period depends on the 
parameters chosen. During the first few Rabi cycles, the 
coherent state parameters a^{t) ~ ag af{t) ~ ae, so 
that the internal and external degrees of freedom appear 
to be decoupled and the system simply oscillates rapidly 
between internal states, this was shown in the first en- 
velope of figure 1. However, for a long time scales, the 
coupling between the internal and external states make 
effects, this result in a modulation of the Rabi oscilla- 
tions. Up to now, the engineered reservoir does not take 
place and the results illustrated in figure 1 indicate that 
without engineered reservoir there are not any decoher- 
ences occur in the system. 

A engineered reservoir coupled to the trapped ion is 
simulated by variations in the trap frequency, oscillating 
near the ion's original trapped frequency. Physically, de- 
coherence in this case arises from random perturbations 
of the Hamiltonian. In what follows, we want to model 
the effects of the variations in the trap frequency with 
the same formalism in ref. [12]. Our main idea is to treat 
the variations in the trap frequency as fluctuations. For 
an ion in a Paul trap [6], Colorado group realizes the vari- 
ations by a random voltage noise source applied to the 
trap electrodes, the noise source is passed through a low- 
pass filter network with a cut-off frequency well below the 
trap frequency. The atom then sees a harmonic poten- 
tial with fluctuating spring constant. The Hamiltonian 
for a trapped ion in a harmonic potential with fluctuating 
spring constant in the rotating frame is 

^i = ^ + Vix)+p^kL^-^^c7, + -mio^e{t)x^, (13) 

this Hamiltonian is just the eq.(7) plus a term which de- 
scribes the fluctuations in the trap frequency. If we take 
the fluctuations as a white noise, i.e.. 



s{t)dt = VfdW{t), 



(14) 



and set 



X 



X, Px = {2mu!) 



H = Lo{P^ + X^) + gP,a, + VfijX^dW{t) 

= Ho + VfojX^dW{t). (15) 

Here, dW{t) is the increment of a real Wiener process[13], 
g = \/2mujk]^.j.^ and P scales the fluctuations. For a 
single run with a known behavior of the fluctuations in 
time, we use a stochastic Schrodinger equation in the Ito 
formalism[14] 



dpi{t) = --[Ho,pi\ 



':^[x\pMw{t)-\[x^[x\pA 

(16) 



to describe the time; evolution of the density operator 
pi in the rotating frame. Here, we are not interested in 
the effects of the fluctuation in short time scale, in this 
sense, wc may take a average over the fluctuation to get 
the master equation for the average density operator 



dp" 
~dt 



= -i[Ho,p'']-V[X^[X\p'^] 



(17) 



We want to determine the off-diagonal element of the 
density operator p = {72P"f^2 ™ internal space 

spanned by |e) and where U2 is given by eq.(4). 
To do this, we first derive a system of equations for 



dp"ge 

dt 



'^igPxPge 



-w^[X^[X\p-V 



-^ = -2igPxP^g-'-u;'[X'[X^ptg] 
^ = -ll^^[X'[X^pt^,{i = e,g) 



(18) 



It is easy to show that the off-diagonal element of the 
density operator p can be represented as 



{g\p\e) = aUt)aUt)p''gg{t) + Km^pUt) 
+ \a2{t)\^p%{t)+al{t)a2{t)pl,{t). 



(19) 



For a short time scale, {n\p1j{t)\m) 0, for m ^ n, The 
module R{t) of the off-diagonal element which character- 
izes the decoherence is given by 



(20) 



the Hamiltonian (14) becomes 



R{t) = Mod{Tre{g\p\e)} 

00 

= Mod{J2Hit)a*^it){n\p''ggiO)\n) 

n=Q 

+ [amf{n\plM\n) 
+ \a2it)\^{n\ptgmn) 

+ at(t)a2(t)(n|p:,(0)|n)]e-r-^("^+"+i)*}, (21) 

where ai{t), a2{t) and |n) are the same as above men- 
tioned, Tvp denotes a trace over the external states. As 
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eq.(21) shows, the decoherence rate rw^(n^ + n + 1) de- 
pends on the character of the fluctuations, the trap fre- 
quency and the motional state of the trapped ion. Physi- 
cally, the trap frequency play a role of the coupling of the 
system to the reservoir, so the larger the trap frequency, 
the larger the decoherence rate(decay rate). The fact 
that the decoherence rate depends on the motional state 
of trapped ion was observed in the experiments of ref. [7] . 
Again, we consider the state given by eq.(lO) as a initial 
state, the numerical results of eq.(20) arc illustrated in 
Fig. 2, there are a few Rabi cycles in the beginning. How- 
ever, for a long time scales, the oscillations disappear, it 
is evidence that decoherence occur in the system. 

To sum up, decoherence in a two-level trapped ion 
is studied in this paper. The decoherence is induced 
by coupling the ion to the engineered reservoir, which 
is simulated by random variations in the trap frequen- 
cies. Without this reservoir, the transitions between the 
ionic internal levels manifest modulated Rabi transition, 
whereas the transition was suppressed when the engi- 
neered reservoir take place. The suppressed transitions 
indicate that there is decoherence in the trapped ion 
system. 
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Figure captions: 

Fig.l: The module of the off-diagonal element of the 
density operator is plotted as a function of time with fixed 
trap frequency. The parameters chosen are = «£ = 3, 
w = 27r X 11.3MHz, S = A.OGHz, VLl ^ lO.OKHz. 

Fig. 2: This plot shows the decoherence in trapped ion 
induced by the engineered reservoir. The module of the 
off-diagonal element of the density operator is calculated 
by the master equation. This plot shows the module as 
a function of time. The parameters chosen are the same 
as in figure 1 and F = l.OKHz. 
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